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£h ■ ABSTRACT 

' We report on the results of a multi-wavelength study of the compact object candidate 

1RXS J141256.0+792204 (Calvera). Calvera was observed in the X-rays with XMM- 
Newton /EPIC twice for a total exposure time of ~ 50 ks. The source spectrum is 
thermal and well reproduced by a two component model composed of cither two 
(absorbed) hydrogen atmosphere models, or two blackbodies with temperatures kT\ ~ 
• 55/150 eV, kT2 ~ 80/250 eV, respectively (as measured at infinity). Evidence was 

(~) found for an absorption feature at ~ 0.65 keV while no power-law high-energy tail is 

£NJ . statistically required. Using pn and MOS data we discovered pulsations in the X-ray 

' emission at a period P = 59.2 ms. The detection is highly significant (> lie), and 

unambiguously confirms the neutron star nature of Calvera. The pulse profile is nearly 
■ sinusoidal, with a pulsed fraction of ~ 18%. We looked for the timing signature of 

' Calvera in the Fermi Large Area Telescope (LAT) database and found a significant (~ 

5u) pulsed signal at a period coincident with the X-ray value. The gamma-ray timing 
analysis yielded a tight upper limit on the period derivative, P < 5 x 10~ 18 Q(E ro t < 
10 33 ergs -1 , B < 5 x 10 10 G for magneto- dipolar spin-down). Radio searches at 
1.36 GHz with the 100-m Effelsberg radio telescope yielded negative results, with a 
deep upper limit on the pulsed flux of 0.05 mJy Diffuse, soft (< 1 keV) X-ray emission 
about 13' west of the Calvera position is present both in our pointed observations 
and in archive ROSAT all-sky survey images, but is unlikely associated with the 
X-ray pulsar. Its spectrum is compatible with an old supernova remnant (SNR); no 
evidence for diffuse emission in the radio and optical bands was found. The most likely 
interpretations are that Calvera is either a central compact object escaped from a SNR 
or a mildly recycled pulsar; in both cases the source would be the first ever member 
of the class detected at gamma-ray energies. 

Key words: Gamma-ray: stars - pulsars: general - pulsars: individual: 1RXS 
J141256. 0+792204 (Calvera) - stars: neutron - X-ray: stars 
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1 INTRODUCTION 

Isolated neutron stars (NSs) have been for a long time as- 
sociated with radio pulsars. It was only during the last 
two decades, thanks to multi- wavelength observations from 
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both space and ground-based instruments, that our pic- 
ture of the Galactic NS population changed. These obser- 
vations unveiled the existence of different types of isolated 
NSs which are radio quiet or have radio properties quite 
at variance with those of ordinary pulsars. Examples are 
the soft gamma repeaters and the anomalous X-ray pul- 
sars (the likely m agnetars; e.g. IWoods fc Thompson! 120061 : 
Mercghetti 2008), the central compact objects (CCOs; e.g. 



IPavlov. Sanwal fc TeteJ 12004 |Pe Lucall2008l ) in supernova 
remnants (SNR), the seven di m ROSAT sources (X-ray 
Dim Isolated NSs, XDINSs; e.g. lHaber]H"2007l : lTurollall2009l ) 
and the more recently discovered rotating radio transients 
(RRATs; [McLaughlin et all 120061 ). Observations at X-ray, 
gamma-ray and optical wavelengths have been crucial in 
revealing the rich variety of NS properties and their phe- 
nomenology. 



Only a limited number 20, in comparison with ~2000 
known ordinary radio pulsars) of sources in each of these new 
isolated NS classes has been identified so far. This makes 
establishing links between different groups hazardous and 
no unitary picture of the Galactic NSs has emerged to date. 
The search for new members of each class is therefore of the 
utmost importance, as is the study of objects with properties 
intermediate between those of known groups, which could 
provide the much needed missing links. 

In this respect, it appears very promising that 
two new, apparently radio quiet X-ray sources have 
been recently proposed as isolated NS candidates. 
These are 1 R XS J 14 1256.0+792204 (dubbed C alvera ; 
Hessels et ail 120071; iRutledge. Fox fc Shevchukl |200S| ; 



Shevchuk. Fox fc RutledgeT 2003 ) , which was f ound in the 



ROSAT Bright Source Catalogue (RBSC, IVoges et all 
I1999D , and, at fainter flux levels 2XMM J104608 .7-594306, 
discovered with XMM-Newton jpires et alj|2009l) . 

Calvera, in particular, is a puzzling source and its in- 
terpretation represented a conundrum (see § [2J ■ At variance 
with most known radio quiet isolated NSs, the source is at 
high Galactic latitude, and it exhibits a quite hot, thermal 
spectrum with little absorption. Calvera was first selected 
in the RBSC and identified as a possible isolated NS can- 
didate on the basis of its large X-ray-to-optical flux ratio. 
The source was then observed with the SVre/J-X-Ray Tele- 
scope (XRT) and with Chandra (see § [2] for details of the 
past observations, and references therein). In particular, the 
first, short (~2 ks) Chandra pointing (obs. ID: 8508) pro- 
vided a refined X-ray position and follow-up Gemini observa- 
tions confirmed the lack of optic al counterparts down to q ~ 
26.3 implying Fx /F opt > 9000 (iRutledge. Fox fc Shevchukl 
2008). However, despite the improved spectral information, 
these data were insufficient to discriminate among the dif- 
ferent interpretations for its nature. 

Clearly, better spectral and timing information is cru- 
cial to unambiguously classify the source. Here, we present a 
multi-wavelength study of Calvera, based on two new XMM- 
Newton observations, on the analysis of publicly available 
Fermi-L AT data, and on a new radio observation taken at 
the 100-m Effelsberg telescope. We also used Chandra and 
ROSAT archival data. The paper is organised as follows. 
We first summarise the results from the past observations 
and the proposed scenarios in § [5] Spectral and timing re- 
sults from the XMM-Newton and Fermi-L AT observations 
are presented in § O while the radio observations are sum- 
marised in §3] A search for potential sources of diffuse X-ray 
or radio emission in the proximity of Calvera is presented in 
§ [5] discussion and conclusions follow in § 161 Tl 



2 PAST OBSERVATIONS AND PROPOSED 
SCENARIOS 

The first X-ray observations of Calvera, taken with the 
Swift-XNT, had very limited counting statistics (74 source 
counts in the 0.3-10 keV range ) and only allowed a very 
prelim inary spectral analysis (|Rutledge. Fox fc Shevchukl 
120081 ): the spectrum appeared rather hard, but the spectral 
model was virtually unconstrained. A blackbody (BB) fit 
gave temperature kTsB ~ 215 eV, radius Rbb ~ 7 km (for 
a 10 kpc distance), while a power-law (PL) model yielded a 
photon index V ~ 2.8. The 0.3-10 keV flux was 1.2 x 10~ 12 
and2.5xl0 -13 ergcm _2 s _1 for the two models, respectively. 
A thermal bremsstrahlung and a Raymond-Smith plasma 
model were also tried, providing temperatures of ~ 0.8 and 
~ 1.5 keV. No estimate of the interstellar absorption was 
derived (in all fits N h was held fixed to th e total Galactic 
value obtained from iKalberla et ail (|2005h ). so the source 
distance was unconstrained. 

On the basis of the very li mited spectral information 
obtain ed with the Swift data, IRutledge. Fox fc Shevchukl 
(2008) discussed several intriguing possibilities for the na- 
ture of Calvera. However, they also concluded that none 
of them was without problems. For instance, they proposed 
that the source may be an XDINS, but hotter than the other 
members of the class (kTsB ~ 40-100 eV). In this case, un- 
der the arbitrary assumption that the emitting radius is the 
same for all objects of the class, the observed flux would 
imply a distance of ~ 2 — 8 kpc from the Sun and a height 
of ~ 1.3 — 5.1 kpc above the Galactic plane. Taken at face 
value, this appears hardly compatible with the star's age 
as derived from its temperature (assuming standard cooling 
models), unless either Calvera has an unprecedented high 
velocity (a few thousand kms -1 ), or did not cool according 
to conventional models, or it has been re-heated by some 
unknown mechanisms. 

A magnetar identification was also recognised to be 
problematic. For an average quiescent magnetar luminos- 
ity of ?a 10 35 ergs -1 , the flux inferred from the best fit of 
the Swift-XRT spectrum would imply a distance of at least 
66 kpc, and a height of 40 kpc above the Galactic plane. 
Calvera would then be the first magnetar to belong to a 
halo population. Although evidence of massiv e compact ob- 
ject f ormation in the halo has been found (e.g. lMirabel et ail 
120011 ). this interpretation was considered statistically un- 
likely. On the other hand, a Galactic plane population ori- 
gin would require a luminosity fainter by about 5 orders of 
magnitude, in which case Calvera would be a candidate for 
the faintest, and maybe the oldest, identified magnetar. 

IRutledge. Fox fc Shevchukl (|2008T ) also considered the 
idea that Calvera might be a CCO. In this case, assuming 
an averag e NS velocity of 400 km s -1 and a cooling age of 
0.5 Myrs (|Page et alj|2004l ). the present location of Calvera 
would imply a distance of ~ 300 pc which would make it a 
factor ~ 10 fainter than other known CCOs. 

The last alte rnativ e proposed by 

iRutledgeT Fox fc Shevchukl (|2008l ). was an old millisec- 
ond pulsar. This scenario would place the source at a 
distance of a few hundreds pc, the exact value depending 
on the assumed size of the hot polar caps. Interestingly, 
for such a small distance, the derived upper limit in the 
optical of g > 26.3 reported by IRutledge. Fox fc Shevchukl 
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(2008) is deep enough to rule out both a white or a red 
dwarf companion. If this is the case, Calvera would thus be 
one of the few solitary field millisecond pulsars discovered 
so far (about 18 known isolated millisecond pulsars in the 
Galactic plane, comprising roughly 30% of the population). 

In order to shed furt h er lig ht on these possibilities, 
IShevchuk. Fox fc Rutled"gel l)2009h re-observed the source 
with Chandra for 30 ks on 2008 April 8 (obs. ID: 
9141). This longer observation allowed them to exclude 
(absorbed) single-component spectral models, as a BB, 
a PL or a pure hydrogen a tmospheric model (NSA; 
IZavlin. Pavlov fc Shibanovlll996l ). since they all gave unac- 
ceptably large values of \ 2 and, in some cases, a value of the 
column density larger than the G alactic one in the so urce 
direction (N gat = 2.65 x 10 20 cm 2 : lKalberla et al-llioolh . In- 
stead, these authors proposed as best fit a NSA model (with 
effective temperature of ~ 110 — 120 eV) combined with a 
spectral feature, either an edge or more likely an emission 
line, at ~ 0.5 — 0.6 keV. The source did not show signs of 
variability neither on short nor long (> 1 year) timescales, 
and pulsations have not been detected down to a 3<r limit 
of 20-30% in pulsed fraction (defined as the semi-amplitude 
of the sinusoidal modulation divided by the mean source 
count rate) in the period range between 0.9 s and 10 4 s. The 
absence of a non-thermal spectral component and the rela- 
tively low limit on the pulsed fraction for periods of a few 
seconds made a magnetar interpretation unlikely, but the 
other possibilities still remained open. Therefore, despite the 
better spectrum, the Chandra data (obs. ID: 9141) did not 
provide sufficient clues to s o lve th e conundrum proposed by 
iRutledge. Fox fc Shevchukl <|2008T ) about the nature of the 
source. 



3 X-RAY AND GAMMA-RAY OBSERVATIONS 

3.1 Source position and association 

We started our analysis by first accurately re-assessing the 
source position determination. This is particularly impor- 
tant because a large positional error precludes firmly ruling 
ou t the association with the close-by (~ 0.9") g ~ 24.8 Star 
A (jRutledge. Fox fc ShevchuklbOOsI ). 

Ba sed on a Chandra /HRC-I observ ation (obs. ID: 
8508), IRutledge. Fox fc Shevchukl (|2008l ) gave: a = 
14 h 12 m 55 i :885, S = +79°22'04'.'10, with a 90% uncertainty 
of 0'.'57, after applying a bore sight correction of Aa — 
-0704 ± 0'.'22 and AS = +0'.'95 ± / .'22 to the Chandra as- 
trometryQ To this aim, they used as a reference the optical 
coordinates of a second source, CXOU J141259.43+791958, 
detected in the HRC-I field, derived from those of its pu- 
tative counterpart identified in their G emini images, cal i- 
brated with the USNO-B1.0 catalogue (|Monet et alll2003l) . 
However, this procedure comes with some caveats. 

T o verify the value of the IRutledge. Fox fc Shevchukl 
(2008) coordinates, we retrieved the HRC-I data set (obs. 
ID: 8508) from the public Chandra archive. First of all, we 
found that their reference X-ray source has a count rate of 

1 As the authors noted, this large correction was probably due 
to a systematic 0'.'4 offset in the absolute Chandra astrometry 
affecting data taken after 2006 December 1. 



« 0.005 cts/s (~ 2.5<r detection significance), which low- 
ers the accuracy of the bore sight correction. Furthermore, 
using only one reference source can introduce systematics 
due to a proper motion of its putative optical counterpart. 
Since the default reference epoch for the USNO-B1.0 coor- 
dinates is 2000. o3 while the epoch of the Chandra /HRC- 

1 observation (obs. ID: 8508) is 2007.13, the effect on the 
bore sight correction can be significant. We identified the 
putative counterpart of the reference X-ray source with the 
USNO-B1.0 star 1693-0051235 (B = 20.51; R = 19.4). Ac- 
cording to the catalogue, this star has a proper motion 
(i a — —6 ± 4 mas yr -1 and (is — 6 ± 1 mas yr -1 , with 
±4 mas yr - bein g a more real istic uncertainty for the lat- 
ter value (see e.g. iGouldl 120031 ). This would imply an off- 
set of Aa ~ +42 ± 28 mas and AS = -42 ± 28 mas in 
the bore sight correction, a ppare ntly not accounted for in 
IRutledge. Fox fc Shevchukl (|2008l ). While this offset would 
be negligible within the error budget of the Chandra and of 
the USNO- Bl.O-based astrometry calibration of the Gem- 
mi images l|Rutledge. Fox fc Shevchukl 120081) , we can not 
rule out that the star has a larger proper motion. We 
queried other astrometric catalogues to verify the USNO- 
B1.0 proper motion. Unfortun ately, the star is too f aint to 
be in the UCAC-3 catalog ue (IZacharias et al.ll2010l) while 
it ap pears in the PPMXL (|Roeser. Demleitner fc Schilbachl 
|2010|) which gives (i a = —7.1 ± 5.6 mas yr and (is = 
6.1 ±5.6 mas yr -1 , consistent with the USNO-B1.0 value al- 
though not statistically significant yet. From this value, we 
set a 3ct upper limit on the star proper motion of (i a = 24 
mas yr -1 and (is = 23 mas yr -1 . Thus, we estimate an addi- 
tional uncertainty of ~ Q' /17 (per coordinate) to the one-sta r 
bore sight correction of IRutledge. Fox fc Shevchukl (2008), 
rising its total uncertainty to 0'.'28 (per coordinate). 

We independently measured the coordinates of Calvera 
in the HRC-I dataset, after it was corrected for the reported 
offset in the absolute Chandra astrometry. Our best fit coor- 
dinates are then a = 14 h 12 m 55!64, S = +79° 22' 03'.' 7, with 
a 90% confidence error of 0'.'6. Given the re-assessed un- 
certainty on the one-star bore sight correction (see above), 
we decided not to apply it to the measured coordinates of 
Calvera. 

A second Chandra observation of Calvera 
was performed with the ACIS-S detector 
l|Shevchuk. Fox fc Rutledgel 120091 . obs. ID: 9141). We 
retrieved the ACIS-S dataset from the public Chandra 
archive and we measured the position of Calvera as 
a = 14 h i2 m 55!84, S = +79°22'03'.'7, with a 90% confidence 
error of 0'.'6. This position is consistent with that measured 
by us using the HRC-I data s et and, obviously, virtually 
identi cal to that obtained by IShevchuk. Fox fc Rutledgel 
(2009) from the uncorrected Chandr a astr ometry. We 
noticed that IShevchuk. Fox fc Rutledgel |2009l) applied the 
bore sight correction of the ACIS-S image, using the posi- 
tion of the optical counterparts to two field sources, CXOU 

2 USNO-B1.0 coordinates are extrapolated to epoch 2000.0 from 
the mean epoch of observations according to the proper motions 
computed by the USNO-B1.0 pipeline. However, the uncertainty 
due to the proper motion extrapolations is not accounted for in 
the nominal coordinate accuracy. This implies larger errors on any 
USNO-Bl.O-based astrometric solution and on the coordinates 
thereby determined. 



4 



CXO/ACIS-S Obs.lD 9141 



. . . r 




MOS 0.2-4.5 keV 



14:13:00.0 14:12:48.0 14:12:36.0 14:12:24.0 18:00.0 16:00.0 14:00.0 12:00.0 14:10:00.0 08:00.0 

Right ascension RA (J2000) 

Figure 1. Left panel: Chandra image of the field taken in 2008 (obs. ID: 9141). Red circles mark the two detected sources: Calvera 
(a = 14 h 12 m 55?83, <5 = 79°22'3'/7, uncertainty = 0.63"at 90% confidence level, cts/s = 0.180 ± 0.003) and CXOU J141252. 41+792152.60 
( a = 14 h 12 m 52?41, <5 = 79°21 m 52'/6, uncertainty = 1.41"at 90% confidence level, cts/s = 0.00052 ± 0.00018). Note that these two 
sources are not resolved in the XMM-Newton data, but, given the low count rate of the second one, we do not expect a significant 
contamination in our spectral and timing analysis. The green circles mark the positions of the two sources used by Shevchuk et al. (2009) 
for the bore sight correction; note they are undetected in our analysis of the Chandra data. Right panel: Combined EPIC MOS1/2 
image in the 0.2-4.5 keV band, with the detected sources marked by red circles. MOS1 was in imaging mode, covering the target and the 
surrounding field (central CCD + 5 outer CCDs), while MOS2 was in timing mode (imaging for outer 6 CCDs). Green circles correspond 
to the green circles on the Chandra image. 



J141220.78+792251.6 and CXOU J141246. 23+792222. 3. 
They then obtained a = 14 h 12 m 55 ! :76, 5 = +79°22'03'.'4 
with a quoted 90% uncertainty of 0'.'31 and 0'.'3 in right 
ascension and declination, respectively. However, we could 
not detect these two sources in the ACIS-S (above a 
detection threshold of 3<r, see Fig. [TJ left panel). We 
confirm, instead, the marginal detection of a faint source, 
CXOU J141252. 41+792152. 60, in the ACIS-S data located 
~ 20" south west of Calvera, which was not detected in the 
HRC-I data (obs. ID: 8508), while the HRC-I source CXOU 
J141259.43+791958 falls out of the ACIS-S field of view. 
Because of the non-detection of t he two sources claimed by 
IShevchuk. Fox fc Rutledgei (2009) in the ACIS-S data, we 
used the nominal Chandra position in the present analysis. 

Finally, for the sake of completeness we checked 
our XMM-Newton observations. XMM-Newton observed 
Calvera twice in 2009, on August 31 and October 10, for 
a total exposure of about 47 ks. The observation s were per- 
formed with the EP IC pn jStriider et al.ll200ll ) and MOS 
jTurner et al l l200lf ) cameras in different CCD read out 
modes (see Table [T]). The net exposure from the first ob- 
servation (obs. ID: 0601180101) was short: less sources were 
detected and no bore sight correction was performed. From 
the second observation (obs. ID: 0601180201), we obtained 
a = 14 h 12 m 55 ! :48, S = +79°22'03'.'3, as computed by the 
EPIC pipeline (lcr statistical error of 0'.'21) which com- 
putes the bore sight correction from the positions of X- 
ray sources in the field matched against those of optical 
sources in, e.g. the USNO-B1.0, NED, and NOMAD cat- 
alogues. In this case, the pipeline finds 35 X-ray sources 
with optical matches. However, we visually inspected the 
X-ray source position on the DSS-2 images and we found 



Table 1. XMM-Newton EPIC observations of Calvera. 



Observation 


EPICM 


Start - End time (UT) 


Exp/ 6 ) 


ID 






(ks) 


0601180101 


pn SW 


2009-08-31 07:14-15:09 


14.0/13.94 




Ml FF 


07:08-15:09 


19.6/19.65 




M2 TU 


07:08-15:05 


20.0/ - 


0601180201 


pn SW 


2009-10-10 04:15-12:27 


14.5/19.48 




Ml FF 


04:09-12:26 


20.3/27.44 




M2 TU 


04:09-12:22 


20.4/ - 



( a ) Instrument configurations: pn SW: Small Window CCD read- 
out mode with 6 ms frame time; MOS1 (Ml) FF: Full Frame, 
2.6 s; MOS2 (M2) TU: Timing Uncompressed, 1.75 ms; The thin 
optical blocking filter was used for all cameras. 
W Net exposures used for timing and spectral analysis (dead-time 
corrected). 

possible matches with 14 optic al sources, of whic h only 4 
are listed bo th in the GSC-2 (fLasker et all l200Sl ) and in 
the 2MASS (|Skrutskie et al.l I2006T ) catalogues. By match- 
ing the pixel-to-sky coordinates of the X-ray sources and 
their associated 2MASS/GSC-2 counterparts we derived an 
astrometric solution with an rms of 0'.'88 and we deter- 
mined a = 14 h 12 m 5S?68, S = +79°22'04'.'2 for Calvera, for- 
mally more consistent with the coordinates obtained from 
the Chandra data sets. 

We notice that t he t wo sources used by 
Shevchuk, Fox fc Rutledgei {2009) for the Chandra bore 
sight are also not detected by XMM-Newton (see Fig. [T] 
right panel). From the sensitivity map of observation 
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0601180201 (pipeline), the 3<j upper limit for the two 
sources is 0.005 cts/s (in the EPIC total band 0.2-12 keV), 
a factor of 30 lower than the count rate of Calvera. This 
converts to a flux of w 3 x 10 -14 ergcm _2 s _1 , with sizeable 
uncertainties depending on the assumed spectral model. 

In conclusion, through our astrometry analysis we con- 
clude that the most reliable coordinates of Calvera are those 
derived from the Chandra astrometry (a = 14 h 12 m 55'!84, 
S = +79°22'03'.'7, with a 90% confidence error of 0'.'6) i.e. 
without applying uncertain bore sight corrections, and they 
indee d confirm that Calvera is not as sociated with star A 
from iRutledee. Fox fe Shevchukl |2008). Moreover, we note 
that the uncertainty on the Calvera coordinates is small 
enough not to hamper a sensitive search for coherent sig- 
nals in our timing analyses (see § 13. 5[) . 

3.2 XMM-Newton data analysis 

For both, spectral and timing analysis we used the lat- 
est XMM-Newton Science Analysis Software (SAS), ver- 
sion 10.0.1. We checked that the pn event files produced 
by "epchain" were clean of unrecognised time jumps before 
we applied an event time randomisation. 

We extracted source spectra for Calvera from the imag- 
ing mode data of pn and MOS1 using a circular region 
around the target position with 15" radius and background 
spectra from a nearby source-free region of same size. Based 
on its low number of counts (see Fig. [T] and caption), 
CXOU J141252.41+792152.60 is not expected to contam- 
inate Calvera's data. Since the spectral response in MOS 
timing mode is currently not well calibrated, we did not use 
the MOS2 timing data for spectral analysis. For spectra we 
selected single pixel events (corresponding to PATTERN=0) 
for pn and all valid events for MOS (PATTERN=0-12), ex- 
cluding bad CCD pixels and columns (FLAG=0). We used 
XSPEC version 12.5.0 for spectral modelling. 

The same pn and MOS event lists considered for the 
spectral analysis were also used as input for the timing 
analysis (with different background screening criteria, see 
Table [TJ. In addition, in order to obtain better temporal 
resolution, we used also the MOS2 data. MOS2 was oper- 
ated in timing mode; in this mode data from the central 
CCD are collapsed into a one-dimensional row to achieve 
a ~1.75ms time resolution. Therefore, instead of a circu- 
lar extraction region, we used a rectangular box of 20 pixel 
width. The local background was measured in a region of 
the images far from the extracted source itself and tak- 
ing care that no other source falls within the extraction re- 
gions by chance (as we already noticed, the Chandra source 
CXOU J141252.41+792152.60 is unresolved from Calvera, 
see § 13. ip . The arrival times of the photons were corrected 
to the barycenter of the Solar System us ing the JPL DE405 
Solar System ephemeris l|Standishl 2004) and the SAS task 
"barycen" . 



3.3 XMM-Newton spectral analysis: pulse 
averaged spectra 

Spectra were accumulated from the pn and MOS1 event 
lists after removing those time intervals affected by so- 
lar proton flares resulting in an effective exposure time of 



~14ks (pn) and ~19.7ks (MOS1; first observation, obs. ID: 
0601180101), 19.5 ks (pn) and 27.5ks (MOS1; second obser- 
vation, obs. ID: 0601180201). 

We started the analysis considering an absorbed dou- 
ble BB model and fitting the two pn and the two MOS1 
spectra simultaneously, allowing only a renormalisation 
factor to account for cross-calibration uncertainties be- 
tween the detectors and possible time-variability of the 
source. To account for interstellar absorption, here and 
in the following we adopted t he elemental abundances of 
IWilms. Allen fe McCravl (|2000l . model "phabs" in XSPEC). 
We note here that the used metal abundances influence the 
derived equivalent hydrogen column density and can intro- 
duce systematic uncertainties when comparing its value with 
Hi-derived column densities. We did not find flux variations 
greater than 3% between the two observations, which is 
consistent, considering the statistical uncertainties, with a 
steady source flux. We then performed an independent fit of 
the spectra from the first and the second observation (still 
joining pn and MOS data within each epoch). The resulting 
best-fit parameters, from the two epochs, were consistent 
within their statistical errors. 

Since the source showed neither significant flux varia- 
tions nor spectral changes between the two XMM-Newton 
observations, in the rest of the analysis we fitted all pn and 
MOS data together in order to obtain the best signal to noise 
ratio. We tested several single, double and triple component 
spectral models, based on different combinations of BBs, 
PLs, NS atmosphere models (NSA) and absorption edges or 
Gaussian lines. In the NSA model, we used R s tar = 12 km 
and M a tar = lAM sun as input parameters for the NS radius 
and mass. The results are summarised in Table [2j where we 
report only the fluxes inferred from the pn data of observa- 
tion 0601180101. As mentioned above, we do not see signif- 
icant flux variations between the two observations: comput- 
ing an average pn flux for the two observations, weighted 
with the exposures, gives a value 0.98 times smaller than 
that reported in Table [2] Also, the fluxes obtained from the 
MOS1 spectra are consistent with those listed in the table. 
We warn the reader that the parameters reported in Ta- 
ble [2] have a different physical meaning inasmuch the best 
fit temperatures obtained by the BB fit are measured at 
infinity, while those from the NSA model are measured at 
the neutron star starface. The two quantities are related b y 
Too = T/(l + z) where l + z = - 2.952M 3tar /i? st(lr is 

the gravitational redshift factor. 

In agreement with IShevchuk. Fox fc Rutlecjge] (2009), 
we find that single component models do not provide a 
good representation of the Calvera spectrum. However, the 
advantage of the EPIC spectra is the higher efficiency of 
the instruments below 500 eV, where the Chandra spec- 
trum has very few counts. Fitting the XMM-Newton data 
with one component models leaves a low energy excess, 
that can be accounted for by a second emission compo- 
nent. Two-component models including a PL (as in the 
combinations PL+BB, PL+NSA) demonstrate a problem: 
the additional PL is always very steep and represents the 
soft part of the spectrum (requiring a high interstellar ab- 
sorption), instead of a hard tail as usually seen from iso- 
lated NSs with a non-thermal component l)De Luca et alJ 
l2005f ). Instead, we found that the best two-component fits 
are composed of either, two thermal components (BB+BB, 
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Table 2. Summary of spectral fits. 



Model 



Xl I d.o.f. N H 



10 20 
cm 




1 843/326 


11.9 


PL + BB 


1 055/324 


8.9±1.0 


PL + NSA 


1.051/324 


9.1+0.8 


BBW 


2.258/326 


0.0 


2BB 


1.053/324 


4.6+1.0 


2BB 


1.096/325 


2 . 7 (/) 


BB X edge 


1.486/324 


0.5±0.3 


BB — gauss 


1.631/324 


0.5±0.3 


2BB X edge 


1.051/323 


2.7(f) 


2BB - gauss 


1.015/322 


2 . 7 C/) 


NSA 


1.347/326 


1.8±0.4 


2NSA 


1.042/324 


6.7±1.3 


2NSA 


1.156/325 


2.7C/) 


NSA x edge 


1.079/324 


2.6+0.5 


NSA + gauss 


1.304/325 


1.5±0.4 


2NSA x edge 


1.053/323 


2.7(f) 



feT Ca) E edge /E line r/a 

eV keV 



7(c) 



erg cm 



06s 

10- 13 

2„-l 



erg cm 



10- 12 

2„-l 



Comment 



3.4 
3.8+0.3 
4.3+0.5 



227+8 
111+7 

198 
93+8/242+9 
110±6/253±10 
201+3 
196±3 
123±11/254±15 
140±14/283±30 

98±3 
27±1/118±I 
671^/150^^ 

97±3 

101±* 
82^ 5 7 /154lf ° 



0.66^17 
0.76±0.02 
0.63±0.03 
0.66±0.05 



0.65±0.02 

0.53") 

0.64±0.03 



0.67+0.08 
0.1(f) 
0.26±0.11 
0.13±0.05 



0.40+0.08 

0.0(f) 

0.31±0.09 



9.12 
8.58 
8.49 

8.12 
8.40 
8.48 
8.30 
8.23 
8.47 
8.51 

8.37 
8.47 
8.56 
8.43 
8.40 
8.53 



0.51 
1.28 

0.84 
0.79/0.67 
0.58/0.59 

1.02 

0.97 
0.61/0.58 
0.82/0.41 

1.65 
1.38/3.00 
0.76/1.16 

1.92 

1.59 
0.46/1.45 



N H > N. 



(e) 
gal 



N H > N gal 
N H > N gal 

N H > N ga i 

EW= -83 eV 
EW= -79 eV 

N H > Ng a l 
EW= 12 eV 



( a ) BB and NSA temperatures are evaluated at infinity and at the NS surface, respectively. It is Too = T/(l + 2), see text for details. 

( b ) Observed flux in the 0.2-10 keV band derived from the EPIC pn spectrum of the first observation. ( c ' Bolometric flux for the 
thermal components (NSA and BB) derived from the E PIC pn spectrum of o bservation 0601180101. (^ No error estimate due to bad 
fit. < e > Ngai = 2.65 X 10 20 cm 2 . lKalberla et al.l ||2005|) . (f> Parameter fixed in the fit. 
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Figure 2. Left panel: EPIC pn and MOS1 spectra fitted with a double blackbody model. Right panel: EPIC pn and MOS1 spectra 
modelled with a double NSA model. Both figures show best-fitting models in which the value of iVjj has been left free. See Table[2]and 
text for details. 



NSA+NSA), or an NSA+edge model. One difference is that, 
while the NSA+edge fit gives a value of Nh consistent 
with the Galactic one (as in the case of the Chandra data, 
IShevchuk. Fox fc Rutled"ge1l2009h , the fits with two thermal 
components give a larger value, which may suggest the pres- 
ence of a local absorption component. However, fixing the 
Nh at the Galactic value does not deteriorate the fits con- 
siderably (see, again, Table[2|. Data and best fitting models 
based on two thermal components (BB+BB and NSA+NSA 
with free Nh) are shown in Fig. [21 Unfortunately, based on 
present spectral data only, it is not possible to discriminate 
between a picture in which the spectrum originates from 
two zones of the NS surface (with temperatures at infinity 



of ~ 80 and 250 eV if the emission is modelled with black- 
bodies, or ~ 55 and 150 eV if NSA models are used), and 
one in which the star surface is at uniform temperature (of 
Too ~80 eV, NSA) with an absorption edge at ~0.65 keV 
present in the spectrum (for a discussion about the possible 
discrimination based on other arguments, see § 16. 1|) . 

We notice that, in the BB+BB and NSA+NSA mod- 
els, the contributions of the two thermal components cross 
near ~0.65 keV, which is probably the reason why the spec- 
tral fit can be accomodated by i ntroducing a feature around 
this en ergy. This is at odds with lShevchuk. Fox &: Rutledgd 
(2009), who proposed as a best-fitting model a NSA plus 
emission line: including a Gaussian emission line in the NSA 
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mo del (with an energy fixed at 0.53 keV, the value derived 
by IShevchuk. Fox fc Rutledgd l2009h improves the fit only 
slightly (see Table [2} and the line equivalent width (EW) of 
12 eV is smaller than that derived from the Chandra spec- 
trum (28 eV). Instead, we found that the fit improves by 
including an additional feature in absorption (the energy 
resolution is not sufficient to differentiate between an ab- 
sorption edge or a negative Gaussian profile). Furthermore, 
we also found that both models based on two thermal com- 
ponents seem to further improve by introducing an absorp- 
tion feature (either, edge or Gaussian line) at ~ 0.65 keV. 
For instance, by adding an absorption edge to the BB+BB 
model (with fixed Nh) gives an F-test statistic value of 7.94, 
and chance probability of 4 x 1O~ 4 

Finally, we tested a 2BB+PL model. The reason for this 
is that there are some data points in the pn spectrum above 
3 keV, although they have low significance (the background 
dominates at these energies, and they could also be caused 
by improper background subtraction). The PL photon index 
was fixed at 2. We found that the PL contribution to the 
total flux is ~ 8%. Again, the fit results in an increased Nh 
as compared to the 2BB case. We conclude that the existence 
of a PL non-thermal component up to a flux level of about 
10% can not be excluded. In order to have further insights 
on this issue, spectra with better signal to noise ratio (S/N) 
above 3 keV are needed. 



3.4 XMM-Newton timing analysis 

Given the rather different observing modes and sampling 
time among the three instruments (At of ~5.67ms, 1.75 ms 
and 2.6 s for pn, MOS2 and MOS1, respectively) we first 
used all the event lists extracted from the three instru- 
ments during the two observations together, and assumed 
the largest (MOS1) sampling time in order to look for sig- 
nals with periods larger than 5.2 s. Moreover, in order to 
maximise the sensitivity of the coherent pulsation search 
we minimise to two (one for each observation) the num- 
ber of averaged power spectra. Significant power spectrum 
pe aks were searc h ed fo r by using the algorithm described 
in llsrael fc Stella! ifl9 96l. No significant peak (above a 3.5<r 
confidence threshold on 2.097.152 frequency trials) has been 
found. The corresponding 3<r upper limit on the pulsed frac- 
tion (defined as the semi-amplitude of modulation divided 
by the mean source count rate) is about 10% in the 5.2-10 3 s 
ra nge (consistent, although tighter, with the constraints set 
bv lShevchuk. Fox fc Rutledgefeoogh . 

The signal search was then carried out by using the pn 
and MOS2 event lists only and with a time resolution of 
5.67 ms in order to sample periods as short as ~12ms. A 
highly significant peak (~11.5cr; see Fig.[3| was found at 
a frequency of 16.89242(2) Hz, corresponding to a period of 
59.19816(7) ms (uncertainties refer to the intrinsic Fourier 
resolution of the power spectrum in Figure[3|. Such a short 




3 We warn that, although widely used, the F-test is not strictly 
speaking a suitable method to test the significativity of (narrow) 
spectral lines and edges. This numbers must therefore be treated 
as an indication only. Given the uncertainties in modelling the 
continuum, further assessment of the significance of the feature 
(e.g. through Monte Carlo simulations) is not warranted. 



0.01 0.1 1 
Frequency (Hz) 

Figure 3. EPIC pn and MOS2 power spectrum together with 
the 4<r detection threshold (dashed line). The power peak corre- 
sponding to the ~59 ms signal is evident. 



and strictly coherent period could be only accounted for by 
the spin period of a rotating NS. 

A more accurate determination of the period, for each 
XMM observation, can be obtained by fitting the phases of 
the modulation over different time intervals. However, we 
noticed that the simultaneous MOS2 and pn light curves 
folded to the above period show a phase shift of about 0.10 
and 0.15 (for the first and second observation, respectively), 
making any results inferred from the phase fitting analy- 
sis unreliable. The presence of both phase shift and pulse 
distortion between MOS and pn was also detected in other 
relatively fast pu lsars, such as PSR B1 706-44 (with a pulse 
period of 102 ms; iMcGowan et al. I l2004h . In order to further 
check for the presence of the observed phase shift we also 
analysed the data from another pulsar, namely PSR B1509- 
58 (pulse period of ~150ms), observed by XMM-Newton in 
September 2000 with the pn and MOS2 in the same obser- 
vational modes also set for our observations. The pn and 
MOS2 folded light curves of PSR B1509-58 are reported in 
Figure!?] The phase shift is statistically significant and equal 
to 0.050±0.001, corresponding to 7.5±0.2ms. Therefore, in 
the following we decide to use only the pn data given that 
no calibration to assess the absolute timing accuracy of the 
MOS exists (see the latest XMM-Newton EPIC technical 
report^fj), and because of the pn larger statistics. 

In order to obtain a refined value of the period we 
divided the first pn observation (obs. ID: 0601180101) in 
four time intervals of duration ~ 5000 s a nd inferred the 
phase of the modulation in each interval (see lDaH'Osso et al.l 
2003;, for more details). The scatter of the phase residu- 
als was consistent with a strictly periodic modulation at 
the best period of P=59. 19821(1) ms or i/=16. 892404(5) Hz 
(90% confidence, 55094 MJD). Unfortunately the inferred 
accuracy is not enough to keep the phase coherency be- 
tween the two observations which are separated by about 
40 days (the inferred period uncertainty of the first XMM- 

4 http:/ /xmm2.esac.esa.int/docs/documents/CAL-TN-0018.pdf 
and http:/ /www2. le.ac.uk/departments/physics/research/src/ 
Missions/ xmm-newton / technical / mallorca-2008-04 / 
08_04_09_mallorca_mos_timode_mgfk.pdf 
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Figure 4. EPIC pn and MOS2 (M2) background-subtracted light 
curves of PSR, B1509-58 folded at the best inferred period of 
~150 ms. The solid lines represents the best fit obtained by adopt- 
ing a model with a sinusoidal plus four harmonics. 
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Figure 5. EPIC pn background-subtracted light curves folded at 
the best periods (see text for details) for the first XMM-Newton 
observation (upper panel, obs. ID: 0601180101), the second ob- 
servation (middle panel, obs. ID: 0601180201) and for the whole 
dataset together (lower panel). 

Newton observation translates to a time uncertainty, at the 
epoch of the second pointing, of ~220 ms corresponding to 
almost four period cycles). Correspondingly, we inferred a 
period separately for the pn dataset of the second obser- 
vation (obs. ID: 0601180201) resulting in a best value of 
P=59. 19822(1) ms or !/=16.892399(4) Hz (90% confidence), 
consistent with the period of the observation 0601180101. 
The corresponding 3<r upper limit on the first period deriva- 
tive P is < 9 x 10~ 15 ss _1 . 

No significant variations have been detected between 
the two observations (see Figure[5)|. The signal shape is al- 
most sinusoidal with pulsed fractions of 18±3% (90% con- 
fidence, defined as (C ma x ~ C m in)/(C m ax + C m i„), where 
Cmax, Cmin are the maximum and minimum number of 
counts). The upper limit on the presence of a second har- 
monic is of the order of 5-7% (3cr confidence level). 

Finally, we searched for pulsation in the Chandra /HRC- 
I archive data (obs. ID: 8508) with the aim of obtaining 



another period measurement and therefore a measure of the 
P or a more constraining upper limit on it. No coherent 
signal was found, and the corresponding 3a upper limit on 
the pulsed fraction is larger than 100%. 

3.5 Gamma-ray data analysis and results 

Gamma-ray observations could represent an important step 
to unveil the nature of Calvera. In particular, the detection 
of pulsed emission at gamma-ray energies can in principle 
be useful to constrain the period derivative (and then the 
rotational energy loss rate) thanks to the large data spans 
typically provided by the current generation of gamma-ray 
telescopes with time tagging accuracy as low as a few /is. 

Hence, we searched for possible gamma-ray counter- 
parts of Ca l vera in the catalogue of gamma-ray sources 
l|Abdo et al.ll2010bl ) detected by the LAT pair- production 
detector on board t he Fermi Gamma-ray Space telescope 
l|Atwood et all [2009). and in the A GILE gamma-ray bright 
source catalogue (jPittori et al.ll2010l ). No entries (spatial de- 
tections) in the gamma-ray catalogues are compatible with 
the X-ray source position within 3cr errors. 

Nevertheless, it was worth performing a timing analysis 
on the LAT data, since the search for source timing sig- 
natures can be more sensitive than the spatial analysis to 
detect and evaluate the flux of weak periodic sou rces (see 
e.g. the case of PSR B1509-58 lAbdo et al.ll2010cT ). though 
typically giving higher statistical errors on the flux. In or- 
der to perform gamma-ray timing analysis we retrieved the 
available LAT public data on Calvera (photon data and 
spacecraft data) through the Fermi Science Support Cen- 
ter (FSSC) web data-server interfacfl We analysed LAT 
data collected from 2008 August 4 (beginning of science 
phase) to 2010 April 27 (LAT runs 239557417 through 
294038606, where the numbers refer to the Mission Elapsed 
Time (MET) in seconds since 00:00 UTC in 1 January 2001). 
During most of this time Fermi was operated in sky scan- 
ning survey mode (viewing direction rocking 35° north and 
south of the zenith on alternate orbits). 

We performed LAT standard data processing using GT- 
SELECT, GTMKTIME and GTBARY tools obtained from the 
HEADAS distribution of the Fermi ScienceTools (version 
v9r15p2) built on Scientific Linux 5 64 bit operating system. 
We selected E>100 MeV events extracted within the "region 
of interest" (ROI) of 2 degrees from X-ray source position 
(S/N ratio of Fermi -LAT pulsars is typically ma ximised tak- 
ing RO I of ~l-2 degrees around pulsar position: lAbdo et al.l 
l2010al ). We selected only events with high probability of be- 
ing photons ("diffuse" event class=3, data quality=l) col- 
lected within 105 degrees maximum zenith angle (applying 
ROI-based zenith angle cut) and excluding data acquired 
during passages through the South Atlantic Anomaly. Solar 
System ephemeris JPL DE405 were used for the barycentric 
corrections. With the above selections, we obtained 2518 
counts for the whole ~1.8 years data span. 

The search for the Calvera timing signature in gamma- 
ray started from the ephemeris provided by the X-ray 
observations (y and upper limit on \v\). Standard epoch 
folding was performed within 10cr errors from the X-ray 

5 http:/ /fermi. gsfc.nasa.gov/ssc/. 
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Figure 6. The v — v contour plot for the the gamma-ray timing 
solutions with Z-test significance >4 sigma. The most significant 
detection is u = 1 6.892401975(2) Hz, u = -1.2(7) X lO" 16 Hz 2 
(timing epoch 55094 MJD). The horizontal dashed line corre- 
sponds to v = 0. See text for details. 



ephemeris. A grid of gamma-ray frequencies and frequency 
derivatives was explored with steps oversampling by a fac- 
tor of 10 the canonical resolution allowed by the data 
span, v res = l/T span = 2 x 10" 8 Hz and u res =2/T^ pan = 
7x 10~ 16 Hz 2 , where T span is the time span of the data. Pear- 
son's x 2 statistics was applied to the 10-bin folded pulse pro- 
file resulting from each period sea rch trial. Bin-indepen dent 
parameter-free Z 2 -test statistics (|Buccheri et al.|[l983l ) was 
also applied to the data. 

A significant pulsed signal from Calvera was detected 
both by x 2 -test and Z 2 -test [y = 16.892401975(2) Hz, 
v = -1.2(7) x 10" 16 Hz 2 , P = 0.059198212396(7) s, P = 
4(2) x 10" 19 fi, Z 2 -test=26.2, cr=4.7 with n=l harmonics). 
A reference epoch of 55094 MJD is assumed for the gamma- 
ray timing analysis. Weighting the corresponding detection 
probabilities with the number of independent v and i> tri- 
als (ritriais ~100) the overall gamma-ray pulse significance 
is ~ 3.7cr, corresponding to a fake detection probability of 
only 3xl0 -4 . We verified that our analysis procedure does 
not produce fake detections even considering v and i> ranges 
much larger (ntriaU > 10 4 ) than those compatible with the 
X-rays ephemeris. In obtaining our timing solutions, the po- 
sition of the source was held at the Chandra coordinates 
reported in § 13.11 We have checked that the positional un- 
certainty (^1 arcsec) does not significantly affect barycen- 
tric corrections and then the rotational parameters resulting 
from our timing analysis. 

Fig. [6] shows the v — v contour plot related to the 
gamma-ray timing solutions with Z-test signficance >4 
sigma. Due to v — i> timing solution degeneracy, as it is 
evident from the plot, positive frequency derivatives cannot 
be in principle excluded although unlikely (we also note that 
the deep optical limit allows to rule out accretion from a bi- 
nary companion) and only upper limits on \v\ < 10 -15 Hz 2 
{\P\ < 5 x 10~ 18 fi) can be conservatively claimed. In case 
of spin-down, this corresponds to a rotational energy loss 
Erot <10 33 ergs" 1 . If the spin-down is caused by magneto- 
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Figure 7. Energy resolved gamma-ray lightcurves (total counts) 
taken with Fermi -L AT. The pulsed signal is detected at ~5cr in 
the whole E>100 MeV band (~6 ms bin resolution). Gamma- 
ray emission from Calvera seems present at least up to ~1 GeV 
though the pulsed signal significance is low at E>500 MeV. 



dipolar braking the upper limit on the magnetic field is 
B <5xl0 10 G. 

Obviously, tight constraints on the ephemeris will rely 
on longer data spans. In particular, >4cr signal-to-noise ratio 
on the i> measurements is expected from 1 year of further 
Fermi-LAT observations. 

The energy-resolved gamma-ray light curves related to 
the most significant solution (corresponding to the cross in 
Fig. [5]) are shown in Fig. [7] A single broad peak is detected in 
the gamma-ray light-curves without significant pulse shape 
variations as a function of energy. 

The pulsed flux was computed considering all the counts 
above the minimum of the light curve, using the expression 
PF = (Ctot - nNmin) / Exp = Cpui/Exp, where C to t is the 
total number of counts, n is the number of bins in the light 
curve, N min are the counts of bin corresponding to the min- 
imum, and Exp is the exposure in cm 2 s units. This method 
is "bin dependent" , but we have checked that different (rea- 
sonable) choices of both the number of bins (i.e., n > 10) and 
the location of the bin centre (10 trial values were explored) 
do not significantly affect the results. The pulsed counts for 
E MOO MeV are C pu i = 455 ± 95 cts (C pui = 275 ± 75 for 
E >300 MeV) corresponding to ~15% of total gamma-ray 
counts (the unpulsed flux can be ascribed to diffuse gamma- 
ray background emission). Comparable results are obtained 
using the the expression adopted for the calculation of the 
X-ray pulsed fraction (C max - C ml „)/(C m ax + C min ). No 
pulsed signal can be significantly (>3cr) disentangled when 
data are restricted to narrower ranges in the low energy 
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part (e.g. 100-300 MeV), although a modulation (in phase 
with the profiles obtained from broader energy bands) is still 
visible in the light curves (Fig. [7). Therefore, only rough 
gamma-ray spectral measurements can be performed with 
the present count statistics. 

We made a LAT exposure cube from the spacecraft 
data file using GTLTCUBE procedure from Fermi Science- 
Tools and created the exposure map using the GTEXPCUBE 
tool. The resulting total exposure of the whole data span 
at the source position for E > 100 MeV is 2.2 x 10 10 cm 2 s 
(5.4 x 10 10 cm 2 s for E > 300 MeV), corresponding to a 
E > 100 MeV pulsed flux PF = (4.1 ±0.9) x 10 -8 ph/cm 2 /s 
(PF = (6.4 ± 1.7) x 10~ 9 ph/cm 2 /s for E > 300 MeV). 
We corrected the pulsed flux accounting for the additional 
source counts falling outside the ROI (2 degrees radius) es- 
timated according to the instrument Point Spread Function 
(PSF). We also verified that the pulsed fraction systematic 
errors due to spectral uncertainties in the exposure calcu- 
lation are below count statistics errors. A comparison of 
E > 100 MeV and E > 300 MeV pulsed fluxes provides 
a rough gamma-ray photon index estimate a = 2.5 ± 0.5 
(100 MeV- 10 GeV). 

In the frame of a spin-powered pulsar interpretation, 
the resulting E > 100 MeV luminosity, assuming a beaming 
angle of 1 steradian, is L 7 = 1.3 x 10 32 d| pc ergs -1 , where 
dkp C is the source distance in kpc. The observed gamma- 
ray luminosity and the upper limit on the spin-down power 
(<10 33 ergs -1 ) constrain the distance below ~ 1 kpc, con- 
sidering a likely gamma-ray conversion efficiency upper limit 
of 109© 



3.6 XMM-Newton Pulse phase spectroscopy 

In order to look for spectral variations in the X-ray flux of 
Calvera we first created background subtracted light curves 
in the energy bands 0.2-0.7 keV and 0.7-2.0 keV and a hard- 
ness ratio by dividing the count rates in the hard band by 
those in the soft band. We used single- and double-pixel 
events (PATTERN=0-4) from the EPIC pn data, combin- 
ing both observations and extracted events from the same 
source and background regions as used for the extraction of 
the spectra. The light curves were folded with a period of 
59.198228 ms, derived from a simple chi-square folding test 
of the data, assuming a constant period. The resulting pulse 
profiles and the hardness ratio as function of pulse period 
are shown in FigH 

We then investigated the possibility of a changing rela- 
tive contribution of the two spectral components with spin 
phase. To this aim, we have produced EPIC pn spectra 
(MOS1 has insufficient time resolution for pulse phase spec- 
troscopy, PPS) around pulse maximum and minimum (each 
covering 0.5 in phase) and fitted the four spectra (two ob- 
servations, two phases) with the double BB model. Simi- 
larly to the case of the pulse average spectra, we fitted all 
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Figure 8. Background subtracted folded light curves obtained 
from the combined EPIC pn data of both observations in two dif- 
ferent energy bands. The pulse profiles are normalised by the cor- 
responding average count rates of 0.197 cts s -1 and 0.133 cts s 
(0.2-0.7 keV and 0.7-2.0 keV, respectively). The bottom panel 
shows the hardness ratio HR = CR{0.7 - 2.0 keV)/OR(0.2 - 
0.7 keV), where CR denotes the count rate in the corresponding 
energy band. 

the pn spectra simultaneously. We performed fits with two 
model flavours: models PPS1 and PPS2. In model PPS1 
only the relative overall normalisations are allowed to vary 
between spectra. This corresponds to a case without varia- 
tion in spectral shape between the two phase intervals. In 
PPS2 we allow different normalisations for all the BB com- 
ponents, i.e. a variation of the relative contribution of the 
two blackbody components is possible. The derived spec- 
tral parameters and bolometric fluxes are summarised in 
Table [3] while spectra are shown in Fig. [9] Again, because 
of no significant long-term flux variations, we only report in 
the table the fluxes inferred from the pn spectra from the 
first observation (obs. ID: 0601180101). 

Both models assume no temperature changes with pulse 
phase, and provide a fully acceptable fit. Given the current 
statistics, we are not in the position to allow more free pa- 
rameters in the fit or to test different spectral models. As one 
can see from Tabled model PPS2 yields a better description 
of the spectral behaviour with pulse phase (F-test statistic 
value of 41.3 and chance probability 6.9x 10 -10 ). This shows 
that the spectral shape changes with pulse phase, as also in- 
dicated by the variations in the hardness ratio (Fig(8]), and 
suggests that the relative contribution of the hotter com- 
ponent is higher during pulse maximum and lower during 
minimum. 



6 An efficiency of ~ 10% is indeed very high, although 
it may be attained in some s ources. As an example, 
Mignani, Pavlov & Kargaltscv (2010) recently revised the dis- 
tance of PSR 1055-52, finding a relatively low value of ~ 300 pc. 
This gives a gamma-ray efficiency as high as 10%, in line with the 
value assumed in this paper. 



4 RADIO OBSERVATIONS 

Following the discovery of pulsations in the X-rays and 
gamma-rays, the X-ray position of Calvera was searched for 
radio pulsations using the Effelsberg telescope with the aim 
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Table 3. Pulse phase spectroscopy 
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Figure 9. EPIC pn spectra from two different pulse phase inter- 
vals covering pulse maximum and minimum from the two XMM- 
Newton observations. The histograms show the double blackbody 
model PPS2 as decsribed in the text and Table [3] 
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of better understanding the nature of the source and in par- 
ticular, in case of detection, to measure the period derivative 
given the higher precision generally achievable in the radio 
band. Observations were taken on MJD 55330 at 1.36 GHz 
with a bandwidth of 240 MHz spread over 410 frequency 
channels using a sample time of 55 [is for 78 minutes and 
for four epochs (a total of 312 minutes). On MJD 55337 a 
further observation was taken using the same setup for 60 
minutes. 

These observations were folded using the X-ray 
ephemeris and searched over a range of frequency dependent 
time delays to correct for the unknown amount of dispersion 
(caused by free electrons along the line-of-sight) between 
Calvera and Earth. These dispersion measures ranged from 
to 1000 cm' 3 pc and the step size was optimised for a 
pulsar with a period of ~ 60 ms. The observations were also 
searched blindly for any periodicities and for single dispersed 
pulses (though impulsive radio frequency interferences at Ef- 
felsberg make this latter type of search very difficult). 

No radio pulses from Calvera were detected down to a 
sensitivity of S = 0.05 mjy for each of the observations as- 
suming a duty cycle of 5%. These observations improved the 
previ ous flux density upper limit (obtained bv lHessels et al.l 
120071 . with the Westerbork Synthesis Radio Telescope) of al- 
most an o rder of magnitude ( by assuming a spectral index 
of 1.7, see Kr amer et al. I [19981) . Assuming a distance d < 1 
kpc, this translates into an upper limit of 0.05 mjy kpc 2 for 
the pseudo- luminosity (defined as S x d 2 ) at 20cm, a limit 
that encompasses 99% of all known pulsars (and 100% of 
the mildly recycled ones to which Calvera could belong, see 



Figure 10. MOS1+MOS2 image in the 0.2-1.0 keV band. The 
spectrum of the extended source west of Calvera has been ex- 
tracted from the region inside the green ellipse. The two green 
dashed circles indicate the regions used to extract the background 
spectrum. Small green and red circles are as in Fig.[T] right panel. 



§ 16.2.20 reported in the version 1.40 of the ATNF Pulsar 
Catalog. 

Note that the Effelsberg observations are only sensitive 
to pulsations and not to continuum e mission; inspection o f 
archival NRAO VLA Sky Survey data ijCondon et alJI 19981 ) 
showed no clear evidence of point-like or extended radio 
sources in the vicinity of Calvera (see § [5j . 
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Figure 11. Spectral fit of MOS1/MOS2 data from the extended 
emission with a NEI model (see text for details) . 

5 SEARCH FOR DIFFUSE EMISSION IN THE 
FIELD OF CALVERA 

In order to search for possible sources of diffuse emission in 
the proximity of Calvera, we analysed in detail the XMM- 
Newton and ROSAT all-sky survey (RASS) images. 

Using XMM-Newton MOS1 and MOS2 data, we cre- 
ated colour images for the second XMM-Newton observa- 
tion. At low energies, different features are clearly visible. 
North of Calvera we can see an area with increased soft 
emission, which however shows a sharp edge at the border 
of the CCD, and therefore it is likely to be related to higher 
noise in one of the MOS CCDs. More interestingly, the im- 
age shows an extended region of soft X-ray emission west 
of Calvera (at the outer rim of the field of view, where the 
source detection software also identifies three sources, see 
the combined MOS 0.2-1.0 keV image shown in Fig. fTU)) . 
The extended emission is relatively soft and best visible be- 
low 1.0 keV. 

The extended emission is also clearly seen in the stan- 
dard RASS images, where it is best visible in the 0.1-0.4 keV 
band, but also seen in the 0.5-2.0 keV band. The ROSAT 
source detection also located a source within the extended 
emission region: this detection is i ncluded in the bright 
source catalo gue dVoges et al.l Il999l ), but lacks an optical 
identification (jzickgraf et al.ll2003l ) The ROSAT images do 
not show significant emission north of Calvera, confirming 
that the feature seen in the EPIC image is likely to be due 
to CCD noise. 

Although it is unlikely that the extended emission is re- 
lated to Calvera, considering their angular separation of 13', 
it remains an interesting detection. Its X-ray morphology is 
reminiscent of an old supernova remnant, or a pulsar wind 
nebula. Further insights on its nature can be obtained by a 
spectral analysis. Therefore, we extracted the MOS1+MOS2 
spectra from an elliptical region (246"and 117" semi-axes) 
covering the extended source. The background was extracted 
using two circular regions (122" radius), free of XMM- 
Newton detections (see Fig. I10[) . Spectra were fitted using 
a non equilibrium ionisation (NEI) model (Fig. Ill[) . which 

7 http:/ /www. atnf.csiro.au/research/pulsar/psrcat/ 



reproduces the X-ray s pectra of faint SNRs. e.g. in the Mag- 
ellanic Clouds (see e.g. lFilipovic et al1l2008l ). We found that 
the spectrum of the diffuse emission is well represented by 
this model (xr = 1-35 for 46 d.o.f.), and clearly shows the 
(unresolved) emission line triplets of O VII and N VI. From 
the fit, O appears to be over-abundant by a factor of ~ 2 — 3. 
Unfortunately, due to the scarce data, the temperature is un- 
constrained (it has been fixed at 0.7 keV during the fit) and 
we could only derive an upper limit on the column density, 
N H < 1.5 x 10 20 cm- 2 (3cr). 

In order to investigate further the nature of the ex- 
tended emission, we inspected the DSS R, B and IR im- 
ages, but did not find evidence for spatial structures. We 
also checked for evidence of diffuse H a emission at the po- 
sition of the XMM-Newton extended source using images 
from the The Virginia Tech Spectral-Line Survey (VTSSfl 
which has pixel size of 1'6 and field of view of 10° x 10° per 
pointing. Unfortunately, the Calvera field is not yet avail- 
able in the processed Survey data. No other H a survey data 
are available, which cover the field with a sufficient spatial 
resolution. 

We then inspected the field in radio continuum using 
data from the NRAO/VLA Sky Survey (NVSS) at a fre- 
quency of 1.4 GHz, which covers the sky at 5 > —40° 
with beams of 45" FWHM (|Condon et al.lll998l ). From the 
NVSS data (Fig. I12[l we could not find evidence for any ra- 
dio source, either point-like or extended, at the position of 
Calvera, down to a flux level of ~ 0.1 mjy (3a) at 1.4 GHz. 
As far as the radio emission from Calvera itself is concerned, 
this upper limit is much less constraining than those dis- 
cussed in §U The closest radio source detected in the NVSS 
data is located at ~ 30" from the position of Calvera, it 
is not extended and has an estimated flux of ~ 0.47 mjy. 
At the same time, we found about a dozen radio sources 
(~ 0.5 — 2.5 mjy) within a circle of 5' radius centred on the 
position of the X-ray extended emission. However, none of 
them shows evidence of extended radio emission on angular 
scales larger than ~ 2' 5. We could not find evidence of dif- 
fuse emission over the whole searched area down to a limit 
of « 1.5 /iJ arcsec -2 , which we assume as the upper limit 
on the radio brightness of the extended emission detected in 
the XMM-Newton and ROSAT data. 



6 DISCUSSION 

We have presented a new multi-wavelength study of the 
NS candidate Calvera (1RXS J141256.0+792204), based on 
two new XMM-Newton observations, on one set of pub- 
licly available Fermi -L AT data and on one new observation 
taken at Effelsberg in the radio band. We also used Chandra 
and ROSAT archival data. We discussed previous determi- 
nations of the source position, and concluded that at present 
the best estimate is that based on absolute Chandra ACIS-S 
astrometry, i.e. a = 14 h 12 m 59!84, S = +79°22'03'.'7 with a 
90% confidence error of 0'.'6. 

The combination of XMM-Newton and Fermi -L AT 
has allowed us to discover the period of the source, P ~ 
59.19 ms, and to place an upper limit on its spin-down 

8 http://www.phys.vt.edu/ halpha/ 
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Figure 12. Left panel: NVSS image of the Calvera field. The contrast has been increased to show faint sources. The position of Calvera 
is marked by the white circle whose radius has been fixed to 10" for a better visualisation. The closest radio source is detected ~ 30" 
northwest of the Calvera position and is well away from the Chandra error circle. The white square (5' X 5') marks the region highlighted 
in the inset (lower left corner). The large white circle (5' radius) is drawn around the position of the extended source seen in the XMM- 
Newton and ROSAT data. No extended radio emission is detected on angular scales larger than ~ 2'5 down to a flux limit of 1.5^tJ 
arcsec -2 . Right panel: Same as in the left panel, with the X-ray contours from XMM-Newton (green) and RASS (cyan) overplotted. 
For XMM-Newton we used the combined MOS 0.2-1.0 keV image (just one contour level is shown, corresponding to 5.4 X 10 — 3 cts s _1 
arcmin -2 ) and for RASS the total 0.1-2.4 keV image (three contour levels: (2.8, 4.1 and 5.5) X 10 -3 cts s — 1 arcmin -2 ). 



rate, \P\ < 5 x KT 18 6. We warn that positive frequency 
derivatives can not be excluded, although, based on the con- 
tours presented in Fig. [B] they appear to be less likely. This 
possibility will be addressed by future timing monitoring, 
and will not be discussed further hereafter. Instead, in case 
of spin-down our measurements correspond to a rotational 
energy loss E rot < 10 33 ergs -1 , a large characteristic age 
> 1.55 Gyrs, and a low magnetic field, B < 5 x 10 10 G 
(under the assumption of magneto-dipolar breaking). 

Although the best-fit of the X-ray spectrum is not 
unique (see Tabled, according to our findings (see § l6.1l for a 
detailed discussion), the most likely spectral model consists 
of two thermal components (blackbodies) with temperatures 
kTi ~ 150 eV, kTi ~ 250 eV, plus possibly an absorption 
edge at ~ 0.65 eV (~ 0.65(1 + z) eV when measured at the 
star surface). Similar spectral features are found in other 
classes of th ermally emitting NS, as the XDINSs (see e.g. 
Turollall2009l. and references therein) or IE 1207.4-5209 (e.g. 
Bignami et ail 120051 ; [Mori fc Hailev||2006l) . which, however, 
are characterised by a much larger spin period. By interpret- 
ing the edge as an electron or a proton cyclotron feature 
gives B = 6 x 10 10 (1 + z) G or B = 6 x 10 13 (1 + z) G, 
respectively: the former value is not too different from 
the limit inferred from the timing (which, however, only 
constrains the large scale dipolar component). Instead, we 
do not find substantial eviden ce for the emission line at 
~ 0.5 — 0.6 keV proposed by IShevchuk. Fox fc Rutledgel 
(2009); on the other hand, considering the thermal char- 
acter of the underlying continuum, a feature in absorption 
appears a more likely outcome. The X-ray luminosity is 
Lx ~ 10 32 d| pc ergs -1 (with small variations depending on 



the assumed spectral model, see Table [2J, and there is no 
evidence for a PL tail up to a flux contribution of ~ 10%. 

Overall, Calvera is a very puzzling pulsar. It has a ther- 
mal soft X-ray spectrum (with a possible absorption fea- 
ture), but, at variance with other similar sources as the 
XDINSs, it is not a "soft" source: through our timing 
search we found that it is still detected at > 4<r even at 
E > 300 MeV. A rough estimate of the 100 MeV-10 GeV 
photon index gives a = 2.5 ± 0.5. The gamma-ray luminos- 
ity is L 7 = 1.3 x 10 32 d| pc ergs -1 , computed by assuming a 
beaming angle of 1 steradian. Assuming a gamma-ray con- 
version efficiency upper limit <10%, and assuming that the 
observed gamma-ray luminosity is powered by spin-down, 
the upper limit on the rotational energy loss (< 10 33 ergs -1 ) 
provides a quite strong constraint on the distance, that turns 
out to be below ~ 1 kpc (much lower distances might result 
from the future estimate of the true value of E ro t). 

The source is not detected in the optical and 
radio bands. Based on the non detection of the 
source on a Gemini-North+GM OS imaging of the field, 
iRutledee, Fox fc Shevch uk (2008) reported a 3<r upper limit 
of g > 26.3 mag, or F g < 0.11 fj.3y, at 4750 A. As for the 
radio band, our Effelsberg observation allowed us to set a 
deep upper limit on the source pulsed emission of 0.05 mjy, 
assuming a duty cycle of 5%. 

Finally, we report evidence for diffuse X-ray emission in 
the field of Calvera, ~ 13' west of the source. 

In the following, we will first discuss the spectral and 
timing results, and then we will elaborate on the implica- 
tions of our new findings regarding the perspectives on the 
source interpretation. 



14 




Figure 13. Left panel: The best-fitting double blackbody model shown in Fig. [2] (left panel) extrapolated in the optical band. Red lines: 
hot component, blue lines: cold component, black lines: total. Solid and dashed-dotted lines correspond to absorbed and unabsorbed 
spectra. The arrow marks the g-band upper limit inferred from Gemini-North+GMOS imaging (see text and Rutledge et al. 2008). Right 
panel: Same as in the left panel for the NSA+NSA model shown in the right panel of Fig. [2] 



6.1 The new timing and multi- wavelength 
spectral properties 

Thanks to the high sensitivity of XMM-Newton at low en- 
ergies, we have been able to reach a good characterisation 
of the X-ray spectrum. In agreement with previous finding s 
based on Chandra data l|Shevchuk. Fox fc Rutledge! i2009). 
we found that the X-ray spectrum can not be reproduced 
by a single component model, while a satisfactory fit can 
be obtained by adding a spectral feature to an atmospheric 
model (NSA). This spectral decomposition has the advan- 
tage that it requires a relatively low value of interstellar 
absorption (Nh ~ 2.6 x 10 20 cm 2 ), compatib le with that in- 
ferred from HI maps in the source direc tion (|Kalberla et al.l 
120051 : IShevchuk. Fox fc Rutledge] 120091) . In this case, the 
star surface emits at a uniform temperature of ~ 100 eV, 
and the feature is more likely to be an absorption edge at 
« 0.65(1 + z) keV. 

The main problem with a scenario based on uniform sur- 
face emission is that our measurement of a relatively large 
X-ray pulsed flux (PF ~ 18%) seems to rule out a relatively 
smooth thermal map. Pulse profiles produced by the ther- 
mal surface distribution induced by a simple core-centered 
dipol a r mag netic field have been investigated long ago by 
IPagel (|l995l ). under the assumption that the surface emits 
(isotropic) blackbody radiation. Because of gravitational ef- 
fects and of the smooth temperature distribution (the tem- 
perature monotonically decreases from the poles to the equa- 
tor), the pulse modulation is quite modest (PF < 10%). 

Indeed, our spectral fits show that the spectral decom- 
position is not unique: an equally likely possibility is that 
the spectrum is made of two thermal components (again 
with possibly a spectral feature at ~ 0.6 — 0.7 keV), re- 
flecting the presence of two zones of the NS surface at 
different temperatures. The two components can be mod- 
elled either with a NSA or a BB, with the difference that, 
since atmospheric models are harder than a blackbody at 
the same temperature, the temperatures inferred from a 
NSA+NSA fit are systematically lower than those obtained 
from a double BB model. As a consequence, a BB+BB and 
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Table 4. Radii of the emitting regions and source distance in- 
ferred from the spectral modelling with two thermal components. 
In the case of the BB+BB fit, the reported values for the ra- 
dius of the hot and cold component correspond to a source dis- 
tance of 1 kpc; for smaller distances they should be scaled down 
accordingly. In the case of the NSA+NSA models, the source 
distance is computed by assuming that the colder component is 
emitted by the whole star surface, and using R B tar = 12 km 
and Mstar = l-4M slln . We warn that, in this case, the inferred 
distances depend on the assumed values, although slightly. For 
instance, would the mass be as large as 2M Bun , the distance will 
turn out to be a factor ~ 15 — 17% smaller. The interstellar ab- 
sorption is left free to vary during the fit or fixed at the Galactic 
value; values are reported in Table [2] 



a NSA+NSA decomposition result in very different predic- 
tions for the amount of flux expected in the optical band 
(see for example Fig. I13|>. and in principle it may be possi- 
ble to use optical observations to discriminate between the 
two options. Unfortunately, in practice, this is unfeasible. 
The g band upper limit from the Gemini-North+GMOS 
imaging translates into ~ 1.5 x 10~ 4 keV cm~ 2 s _1 keV _1 
at ~ 2 x 10~ 3 keV, which is still compatible even with the 
large fluxes expected in the atmospheric emission scenario 
(see Fig. 1131 right panel). In order to use this approach to ro- 
bustly discriminate among X-ray spectral models, we would 
need an upper limit on the optical flux at least 1-2 orders 
of magnitude lower, which requires deep imaging of the field 
down to g ~ 29 — 31. This is clearly unfeasible with the 
current instrumentations. 
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In this respect, some insights may be obtained exploit- 
ing the detection of the source at high energy based on the 
LAT timing analysis. In fact, the fits to X-ray data with two 
thermal components using a NSA+NSA or a BB+BB model 
result in different estimates of the distance and of the size 
of the emitting regions (see Table 2}. In Table HI we report 
the predicted distances inferred by the fits based on a dou- 
ble NSA model, for the limit case of the colder component 
emitted from the whole NS surface. This is an overestimate 
for the area of the colder region, since a second, hotter NSA 
component is present. However, the ratio of the cold/hot ar- 
eas is « 0.02 — 0.05 (similar to that implied by the BB+BB 
fit), and our estimates of the distance would not change sig- 
nificantly by assuming that ~ 95% (instead of 100%) of the 
whole surface is emitting at the lowest temperature. In this 
case, if Nh is fixed at a value compatible with the Galactic 
absorption, the distance turns out to be quite large, ~1.5- 
2 kpc, which seems too large when compared with the limit 
inferred from the gamma-rays (unless the gamma-ray con- 
version efficiency is extremely high). Similarly, the BB+BB 
model is only compatible with emission from the whole sur- 
face (from the coldest component) if the star is at 5 kpc or 
more, which, again, contradicts the gamma-ray limit. 

In fact, a BB+BB model is compatible with the rela- 
tively large PF only for a d ~ 1 kpc, which gives quite small 
BB radii (~ 0.5 — 3 km) suggesting a scenario in which 
the two thermal components originate in two different small 
spots at the star surface. If this is the case, the rather small 
emitting areas, together with a pulsed fraction of ~ 18%, 
and the sinusoidal pulse shape, are suggestive of a geom- 
etry in which the star is a quite aligned rotator seen at a 
large inclination angle. This is also supported by the pulse- 
phase spectroscopy and hardness ratio variation, which in- 
dicate that there is only a moderate spectral evolution with 
phase. This may be explained if both emitting caps are al- 
ways (partially) in view as the star rotates, with the hotter 
spot becoming more visible near the pulse maximum. The 
gamma-ray pulses are also quite sinusoidal, with no evidence 
for changes at different energies and a pulsed fraction similar 
to that in the X-rays. However, the current impossibility of 
phase-align the X- and gamma-ray observations (see § l3.4l for 
all details) prevents to reach any conclusions on the relative 
positions of pulse maxima. 



6.2 The nature of the source 

Our new findings, in particular the new discovery of the X- 
ray and gamma-ray period and the measurement of an upper 
limit for the spin-down rate, unequivocally identify Calvera 
as a relatively fast spinning NS, and provide long-sought 
crucial information to shed light onto the conundrum about 
its nature (see below and § [5J ■ 

Based on the p reliminary spectral informa t ions a vail- 
able at that time, iRutledge. Fox fc Shevchukl (120081 ) at- 
tempted to classify Calvera by comparing it with differ- 
ent classes of compact objects that occupy different re- 
gions in the luminosity /effective temperature/emitting ra- 
dius spaces. In light of the new discovery of a relatively fast 
period, some of the previously proposed scenarios as that 
of an XDINS or a magnetar are now ruled out. Similarly, 
Calvera's spin period is too large to be compatible with an 



15 




10 3 0.01 0.1 1 

P (s) 



Figure 14. P — P diagram for Galactic field NSs: the red star 
marks our upper limit for Calvera, black dots are radio pulsars 
(from the ATNF pulsar catalogue; Manchester et al. 2005), the 
blue squares are CCOs with a measured (or an upper limit value 
of) P (Gotthelf & Halpern, 2007, 2009; Halpern & Gotthelf 2010) 
and violet triangles are magnetars (data from Mereghetti, 2008). 
The NSs encircled in green are sources detected by Fermi (Abdo 
et al. 2010a). Dashed lines indicate equal magnetic field, while 
dotted lines indicate equal value for the rotational energy losses, 

Erot- 

interpretation in terms of a (fully) recycled millisecond pul- 
sar. 

In the context of the P-P diagram, Calvera appears as 
one of the transition objects that populate the zone between 
the bulk of the pulsar crowd and the group of recycled mil- 
lisecond pulsars (in the bottom-left corner of Fig. I14[) . Also, 
the rotational parameters (hence the values of E and B) are 
consistent with those of the CCOs, the most similar of which 
is the source in Kes 79. 

These findings, taken together with the fact that the 
source is apparently isolated and it is not associated with 
a SNR (but see fl6.2.1[) . open exciting novel possibilities, 
among which the most likely are that of an "orphan CCO" 
or that of a mildly recycled source that was once member of 
a high mass binary system (see § 16.2.11 16.2,2[) . 

6.2.1 A new CCO? 

Based on the timing properties and on the thermal char- 
acter of the spectrum, a first possibility is that Calvera is 
a new CCO. In this respect, we note that the X-ray en- 
ergetics of the source is also unusual. The upper limit on 
E rot (< 10 33 ergs -1 ) is only a factor of 10 larger than 
the X-ray luminosity inferred from XMM-Newton data 
(~ 10 32 ergs -1 at 1 kpc). While an efficiency of ~ 10% 
in gamma-ray is not unusual (see £|3.5f) . the conversion ef- 
ficiency in the X-ray band is typically lowe r, of order of 
~ 10~ 2 - 10" 3 ([Becker. Huang fc Prinzll2010h . If the X-ray 
emission if powered by the source spin-down (for instance in 
the case in which it originates from two hot spots at the star 
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surface, heated by back-flowing magneto-spheric currents), a 
similar efficiency would imply a very low distance, of order of 
~ 100 — 300 pc (or even less, considering that our estimate 
of the upper limit on E ro t is conservative). Alternatively, 
the thermal X-ray emission may have a different origin (e.g. 
anisotropic surface cooling). In this case we note that, based 
on our spectral fitting, we cannot exclude that a ~ 10% of 
the X-ray flux is ascribed to a non thermal PL component; 
if this fraction of the total X-ray emission is spin-down pow- 
ered, an efficiency of 10 -2 translates in d ~ 1000 pc. 

Interestingly, apart from binary systems and magne- 
tar sources, that are thought to be powered, respectively, 
by accretion or by the super-strong magnetic field, the 
only sources that may be charac terised by E ro t 5, Lx 
are t he (slowly rotating) XDINSs l)Kaplan fc van Kerkwiikl 
120091 ) and the CCOs. For instance, the only CCO for 
which the spin period derivative has been (recently) 
measured (i.e. PSR J1 852+0040 in Kesteven 79, see 
lHalpern fc Gotthel j [20ich . has L x ~ lQE Tat . In all other 
cases only an upper limit on the period derivative is avail- 
able; still, IE 1207.4 -5209 has L x ~2 x 10 33 ergs _1 and 
E rot < 10 32 ergs- ^Gotthelf fc Halpernl 120071 ). Puppis A 
has E ro t < 10 33 ergs -1 , that does n ot exclude an efficiency 
of~ 10 -2 (|Gotthelf fc Halperij2009h . The X-ray emission of 
CCOs is thermal and characterised by a double black body 
spectrum, without the presence of hard X-ray tail, similar 
to the X-ray spectrum Calvera. 

A possible counter-argument for the CCO interpre- 
tation might be raised b y com paring with Fig. 2-3 by 
iRutledge. Fox fc Shevchukl (J2008) : Calvera's surface tem- 
perature is quite low with respect to those typically observed 
in this class, and, since the typical X-ray luminosity of other 
CCOs is ~ 10 33 ergs -1 , this interpretation would require 
a source distance of 2-3 kpc (against the constraint from 
the gamma-rays). Nevertheless, we warn that these prob- 
lems with a CCO interpretation may be not too severe. The 
comparison of Calvera with CCOs based on the tempera- 
ture comes with caveats, since it does not account for the 
fact that different CCOs might have different ages, or that 
the thermal history of a CCO might have been modified by 
a recent phase of accretion. Also, we note that the derived 
CCO luminosities are based on the estimated distances of 
the host SNRs, which are usually affected by uncertainties 
of the order of 50% or more. Thus, using CCO luminosities 
as standard candles could be risky. 

The possible association with a SNR, or lack thereof, is 
a crucial piece of information to determine whether Calvera 
is, or not, a CCO. We already noticed that the extended 
X-ray source detected west of Calvera in the XMM-Newton 
and ROSAT data shows evidence for triplets of O VII and 
N VI in the spectrum, which are robust signatures of SNRs 
(although a final confirmation could only come from Ha/S 
II/O III ratios). If we assume that Calvera has moved 13' 
away from the centre of this putative SNR, the kinematic 
age of Calvera (and in turn the age of the SNR) results t = 
39000d fcpc /i>ioo yrs, where vioo is the source velocity in units 
of 100 kms -1 . This is much lower than the 1.55 Gyr age 
estimated from the spin parameters. However, this is a com- 
mon characteristic of CCOs: the lower limits on their spin- 
down ages (although smaller than that of Calvera) are al- 
ways much larger than the SNR ages. A way out, within th e 
"anti-magnetar" interpretation (|Hal P ern fc GotthelJlioioh . 



is to assume that the initial spin periods of these pulsars 
were very close to their current values, so that they have not 
changed much during the pulsar's lifetime. In other words, 
the spin-down age is not representative of the true age of 
the source but only reflects its initial conditions at birth. 

Despite these considerations, we consider it unlikely 
that this extended emission is associated to Calvera. At a 
distance of ~ 1 kpc a typical CCO-associated SNR is ex- 
pected to be much larger: for instance the 7 kyr old SNR 
associated to IE 1207.4-5209 would appear with an angular 
size of ~ 2° /dkpc- Instead, the size of the extended emission 
derived from the ROSAT contours is 16' x 8.5', which sug- 
gests that the SNR is a background object. We note that 
the upper limit on the column density inferred from the 
SNR spectrum (N H < 1.5 x 10 20 cm~ 2 ) is smaller than the 
Galactic Nh in the direction of Calvera. This would make 
the extended emission likely closer, however this comparison 
is hampered by the limited photon statistics of the diffuse 
source. Considering the angular separation and the low E ro t 
of the pulsar, it is also unlikely that the extended source is 
a pulsar wind nebula associated to Calvera. 

On the other hand, if Calvera is an underluminous CCO 
at a low distance of only ~ 300 pc, the large spatial scale 
would make it hard to detect the SNR with narrow field 
X-ray instruments like EPIC. Furthermore, it can not be 
ruled out that Calvera is an old 1 Myr) CCO and that 
its host SNR has already expanded and faded away in the 
interstellar medium. 

This "orphan CCO" scenario is potentially intriguing, 
especially in view of the possible connections between old 
CCOs and other NS classes. Still, in this scenario Calvera 
would be a quite extreme case. So far, no CCO has been 
detected at gamma-ray energies. This may be due to obser- 
vational limits: firstly, in cases in which the embedding SNR 
dominates it is difficult to resolve the gamma-ray emission 
from the CCO. Secondly, detecting CCOs as weak gamma- 
ray pulsars is difficult because only three of these sources 
have an X-ray period and only one (Kes 79) has a measured 
P. In Fig. [15] we report the ratio E ro t/d 2 (a often used fig- 
ure of merit for the detectability in the gamma-ray band o f 
a source with E ro t < 10 34 erg s -1 , e.g. lAbdo et al.ll2010ah 
plotted against the spin period for various relevant samples 
of neutron stars (see caption of Fig. I15[) . As we can see, 
among the CCOs with a measured (or an upper limit value 
of) P, two of them have a relatively low value of E ro t/d 2 , 
which may hamper the gamma-ray detection. Would a CCO 
nature be confirmed, Calvera will represent a unique case 
and our timing technique, here applied to Fermi data, con- 
stitute a key tool to investigate the high energy spectrum 
of other members of the class (the most promising candi- 
date for this search being Puppis A, whose upper limit for 
Erot/d 2 is, as in the case of Calvera, compatible with rela- 
tively a large value). 



6.2.2 Is Calvera the first discovered mildly recycled 
gamma-ray pulsar? 

The spin period and the upper limit on the spin period 
derivative place Calvera well within the area of the P — P 
diagram which is known to be populated by the so called 
mildly recycled radio pulsars. Hence the plausible and very 
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Figure 15. Ratio E ro t/d 2 plotted versus spin period for various 
samples of NSs. The green filled dots represent all the pulsars 
detected so far in gamma-rays by Fermi -L AT (according to Abdo 
et al. 2010a). The black empty dots represent all the radio pulsars 
in the ATNF pulsar catalogue having spin period between 20 and 
100 ms and surface magnetic field smaller than 10 11 Gauss (i.e. 
those generally interpreted to be mildly recycled pulsars; none of 
them has been detected so far by Fermi -LAT in gamma-rays). 
The red star represents the position of Calvera, based on our 
upper limit on E ro t and on a (low) distance of 100 pc. The blue 
filled squares refer to the three CCOs with a known value (or an 
upper limit) for the spin period derivative and hence for E ro t- 
The black cross reports the position of the mildly recycled pulsar 
in the Double Pulsar binary PSR J0737-3039, for which there 
is a tentative detection (to be confirmed) of pulsed emission in 
gamma-ray with Agile-GRID (Pcllizzoni, in preparation), and no 
published detection with Fermi-LAT. All the distances in this 
plot (but for the three CCOs) have been derived from the best- 
estimate reported in the ATNF pulsar catalogue at 30 June 2010 
(Manchester et al. 2005). In most cases, they are inferred from 
the dispersion measure of the radio pulsar. 

intriguing possibility that Calvera represents the first case 
ever of a mildly recycled gamma-ray pulsar. 

It is commonly believed that recycled radio pul- 
sars are formed when a relatively old neutron star 
is spun up by accretion while in a binary system 
( Bisn ovatvi-Kogan fc Komberdll974l ). The amount of mass 
accreted, hence of spin up, is related to the duration of the 
X-ray binary phase, and ultimately depends on the mass of 
the companion: low mass X-ray binaries are thought to be 
the progenitors of the so called "fully" recycled radio pulsars 
(often referred to as millisecond radio pulsar^]) with spin pe- 
riods of few ms, while binary systems with an intermediate 

9 By comparing again with Fig. 2-3 in the paper by 
iRutledge, Fox fc Shevchuk! <200d) . we can see that in principle, 
if the source distance were low (~few 100s of pes), the properties 
of Calvera inferred from the BB+BB fit (i.e. typical temperatures 
around 100-200 eV and rescaled blackbody radii of ~ 0.1 — 0.2 km) 
might be compatible with those of the fully recycled radio pul- 
sars observed in 47 Tucanae. However, the now measured spin 



or an high mass companion give rise to mildly recycled radi o 
pulsars with periods of tens of ms (e.g. lLorimer et alj|2004h . 

In this picture we note that the lack of a detection of 
Calvera in the radio band (down to a very deep flux density 
limit) can be plausibly explained by the relatively narrow 
radio beam missing the Earth. As to the X-ray luminosity, 
one can apply similar considerations to those reported for 
the CCO hypothesis (see § I6.2.1[) . i.e. a rotation-powered 
nature requires a distance as low as ~ 100 — 300 pc. The 
double BB X-ray spectrum, however, is somehow unusual 
for a recycled pulsar, since most of them are characterised 
by either single component (either BB or PL) or BB+PL 
spectra. 

The observability in the gamma-ray band is an unprece- 
dented feature for NSs with the rotational parameters of 
Calvera and one may wonder if it is due to some pecu- 
liar emission process operating in this source or to obser- 
vational biases. In this respect we first note that, provided 
the distance is around or below 300 pc, the gamma-ray ef- 
ficiency nicely fits with the average value (around few per- 
cent) seen in the gamma-ray pulsars observed by Fermi- 
LAT. Even assuming a larger distance of 1 kpc, the gamma- 
ray efficiency remains compatible with what observed in a t 
least one other object (|Mignani. Pavlov fc K argaltscv 2010h . 
Also, the gamma-ray spectrum is within the range of those 
shown by the Fermi pulsars, although pointing toward the 
tail of the steepest spectral indices in the spectral parameter 
distribution. By inspecting Fig. \T5\ we can see that the large 
majority of the NSs with rotational parameters in the range 
of those of the mildly recycled radio pulsars (empty dots in 
the diagram) have a E rot /d 2 ratio well below the minimum 
valu e which resulted in a detection with Fermi so far (filled 
dots. lAbdo et al ]|2010al l. Although the location of Calvera in 
this plot is highly unconstrained, we notice that in principle 
its E ro t/d 2 ratio can be at least ~ 10 times larger (for a dis- 
tance less than 300 pc) than that of the best cases among the 
mildly recycled radio pulsars (similar results hold when plot- 

ting t he analogous parameter sfkZtld 2 , e.g. |Abd o et al.l 
l2010al ). Therefore, in the framework of the interpretation of 
Calvera as a mildly recycled gamma-ray pulsar (for which a 
short distance is largely preferred), it seems that the rela- 
tively large observed gamma-ray flux of Calvera is due to its 
proximity rather than to a stronger intrinsic emission with 
respect to the other known rotation powered neutron stars. 

It is worth noting that the relatively steep gamma- 
ray spectral index of Calvera could be echoed in the spec- 
trum of the mildly recycle d pulsar in the Double Pulsar 
binar y (PSR J0737-3039A, iBurgav et al.ll2003l ; ILvne et all 
2004). For the latter (whose E/d 2 is also very promising, see 
Fig. ll5|) there is a tentative detection with Agile-GRID (Pel- 
lizzoni, in preparation) and no detection so far with Fermi- 
LAT, which would indicate a very soft gamma-ray spectrum 
(interestingly enough, also the X-ray spectrum of the Double 
Pulsar appears unusually soft among the recycled pulsars, 
IPellizzoni et alJl200Sl . iPossenti et"afll2008t) . If confirmed by 
future analysis, this may be a peculiarity of the high energy 
emission from the class of mildly recycled pulsars. 

As to the evolution of Calvera, the presence of a main 



pe riod rules out the millise cond radio pulsar scenario proposed 
bv lRutledge, Fox fc Shevchukl J200g|) . 
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sequence, or a giant or a white dwarf companion can be ruled 
out by the deep optical upper limit. This in turn suggests 
that Calvera is not the descendant of an intermediate-mass 
X-ray binary, since the end products of these systems usually 
comprise a neutron star orbiting a heavy Car bon-O xygen 
white dwarf or (in some peculiar cases, see e.g. Il1|2002[ ) a He- 
lium white dwarf. Given the considerations above, it seems 
very likely that Calvera was recycled in a high-mass X-ray 
binary, in which also the companion star eventually experi- 
enced a supernova explosion. The most probable outcome of 
this event is the disruption of the binary, with the release of 
two isolated neutron stars, one having the typical parame- 
ters of a ordinary pulsar and the othe r - usually dubbed d is- 
rupted recycled pulsar, DRP (i.e. see iLorimer et al.ll2004h - 
having a moderate spin rate (in the range of few tens of ms) 
and a surface magnetic field which is intermediate between 
that of the ordinary pulsars and that of the fully recycled 
pulsars. Alternatively, there is the possibility that the binary 
system survives the supernova explosion, leading to the for- 
mation of a double neutron star (DNS) binary. If the absence 
of a binary companion will be supported by future dedicated 
campaigns of observational. Calvera will enter the rare class 
of the DRPs, which would make its discovery particularly 
worthwhile for investigating various still open issues on the 
final stage of the evolution of high-mass X-ray binaries and 
on the supernova kick. 

For instance, there is currently a large mismatch be- 
tween the theoretical expectations on the relative number 
of DRPs and DNSs and the result of the observations. Ac- 
cording to the fact that the survival of the binary in the 
second supernova explosion requires properly tuned param- 
eters for the pre-supernova system and/or for the kick as- 
sociated to the supernova event, population synthesis stud- 
ies (e.g. IPorteeies-Zwart &: Yungelsonl 1 19981 ) indicate that 
DRPs are expected to be generated at a significantly higher 
rate than DNSs. Since the rotational parameters of the NSs 
are similar in the two classes of objects, their lifetime in the 
radio band and their radio emission properties should also 
be similar and hence one wou ld expect to detect u p to ~10 
times more DRPs than DNSs l|Lorimer et al.ll2004h . At vari- 
ance with this prediction, to date 8 DRPs in the Galactic 
field are reported, a number matching that of the 8 known 
DNSs (from ATNF pulsar catalog at 30 June 2010). Super- 
nova kicks smaller than those usually assumed in population 
synthesis have been proposed to at least partially alleviate 
this problem: in fact that would reduce the gap between the 
birth-rates of the DRPs and that of the DNSs. 

Notably, the distances (derived from the dispersion 
measure) of the known DRPs (all selected in the radio band) 
cluster mostly between 2 and 3 kpc, with the closest source 
(PSR J2235+1506) located at 1.2 kpc. Therefore, Calvera 
would turn out to be much closer than all the other DRPs, 
suggesting the existence of a significant (if not even domi- 
nant) contribution of the gamma-ray pulsars to the popu- 
lation of the mildly recycled neutron stars. Were this the 
case, the current observed ratio between the number of ra- 
dio selected DRPs and DNSs could not be representative 
of the whole population. Of course, additional discoveries 



Our present data do not allow us a search for orbital period- 
icities. 



of mildly recycled gamma-ray pulsars (with particular em- 
phasis for blind searches), as well as detailed binary pul- 
sar population synthesis models are necessary to assess the 
biases due to the small-number statistics and to properly 
test the aforementioned hypothesis. The confirmation of a 
DRP interpretation for Calvera and the discovery of similar 
sources would have strong implications: e.g. it will impact 
on the estimates of the beaming factor (i.e. the fraction of 
sky swept by the emission beam(s)) in the radio and gamma 
bands for the mildly recycled neutron stars, thus constrain- 
ing their emission models; it will lead to increase the overall 
birth-rate of the descendants from high-mass X-ray binaries 
in the Galaxy; it will also suggest the existence of some DNS 
binary at close distance from the Earth and which escaped 
radio detection so far. This may in turn lead to an upward 
revision on the expected rate of events of merging of DNSs 
in the Galaxy, a key prediction for the current generation of 
ground-based gravitational wave detectors. 

Once the spin down rate (and hence the surface mag- 
netic field) of Calvera will be measured it might be possible 
to estim ate its initial spin p eriod (post spin-up phase, see 
Fig. 2 in ILorimer et alj |2004). as well as the time since the 
occurrence of the second supernova explosion in the pro- 
genitor binary. With the present data we can only notice 
that, in both the DRP and DNS hyphotheses, Calvera is ex- 
pected to be m uch younger than the age inferred from the 
spin-down rate llArzoumanian. Corder fc Wassermanl Il999l : 
ILorimer et al] |2004D: e.g. for our upper limit on P, the time 
elapsed from the second supernova would be 10 8 yr, under 
the hypothesis that Calvera was spun up by accretion of 
mass at the Ed dington rate up to t he limit imposed by its 
magnetic field, |Lorimer et alj|2004l ). Since the position of 
Calvera - slightly outside the Galactic disk, at a distance 
< O.Gdkpc above the Galactic plane - could be ascribed to 
the kick imparted to the neutron star (or to the DNS binary) 
by the second supernova, a constraint on the time since the 
explosion (even better if complemented with constraints on 
the proper motion of the source) will in turn allow one to 
study the kick imparted to the neutron star (or to the DNS 
binary) in the supernova. 



7 CONCLUSIONS 

Thanks to our multi-wavelength campaign, we have been 
able to recognise in Calvera an intriguing, low magnetised 
pulsar with a relatively fast period, ~ 59 ms. This rules out 
most of the previously proposed scenarios for the nature of 
the source. Calvera's properties are reminiscent of those of 
CCOs and mildly recycled pulsars, although other interpre- 
tations can not be ruled out (mainly because our timing 
solution is still compatible with a spin up scenario). On the 
other hand, even in comparison with the other members of 
these two classes, the source can be singled out thanks to its 
unique characteristics, mainly the hard gamma-ray emission 
joined to a thermal X-ray spectrum. 

The discovery of peculiar, possibly transitional, sources 
such as Calvera is of the utmost importance to gain a full 
understanding of the Galactic NS phenomenology and to 
eventually establish links between the different classes into 
which isolated NSs have been catalogued up to now. The 
existence of bridging objects is progressively coming into 
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view, as in the case of PSR B1509-58 for the magnetar/high 
field spin-powered PSRs connection llPellizzoni et al. 



of the newly discovered radio magnetar (ILevin et al.112 010). 
or th e bursting young pulsar PSR J1846-0258 (|Gavriil et al.l 
2008). Furthermore, although Calvera is most likely to be 
isolated, the current optical limits do not allow us to exclude 
a NS companion, in which case the source would represent 
one of the rare DNS systems discovered in the Galaxy. 

Given the unique placement of Calvera in this frame- 
work, further investigations aimed at better assessing its 
properties are definitely warranted. Accurate X-ray timing 
can provide a positive determination of the source spin-down 
rate, and a better characterisation of its spectrum, especially 
concerning the presence of a high-energy PL tail and of spec- 
tral features, is bound to reveal much on the processes which 
power its X-ray emission, and ultimately on its true nature. 
Presently we were able only to place an upper limit on the 
pulsed radio flux. Calvera may be genuinely radio quiet but 
the search for radio emission must be pursued further since a 
detection at radio wavelengths would yield an independent, 
accurate determination of P, the dispersion measure and 
the star proper motion (the latter may be also p ursed in the 
X-ray band, see lShevchuk. Fox fe Rutledgel2009h . Follow-up 
X-ray and optical narrow-band imaging observations will al- 
low one to better characterise the properties of the putative 
SNR, and to constrain both its age and distance. All these 
issues are crucial to assess if there is indeed a connection 
between Calvera and the nearby diffuse emission, the pres- 
ence of which was first reported in this investigation. New 
observations at gamma-ray energies will allow one to bet- 
ter constrain the source spin evolution, shed light on the 
spectral distribution at gamma-ray energies and how the 
latter changes with rotational phase. Finally, the search for 
gamma-ray emission from other sources with similar prop- 
erties (CCOs, mildly recycled and low E rot PSRs, ...) using 
the novel timing technique which led us to the detection of 
Calvera in the gamma-rays is certainly in order and will be 
matter of future work. 
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